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Constraints on decaying Dark Matter from XMM-Newton 
observations of M31 
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ABSTRACT 

We derive constraints on the parameters of the radiatively decaying Dark Matter (DM) parti- 
cle, using the XMM-Newton EPIC spectra of the Andromeda galaxy (M31). Using the obser- 
vations of the outer (5'- 13') parts of M31, we improve the existing constraints. For the case 
of sterile neutrino DM, combining our constraints with the latest computation of abundances 
of sterile neutrinos in the Dodelson-Widrow (DW) scenario, we obtain the lower mass limit 
nis < 4 keV. which is str onger than the pre vious one m, < 6 keV. obtained recentlv bv Asaka 
et al. (|2007[). Comp aring this limit with the most recent results on Lyman-o- forest analysis 
of IViel et alj (l2007h (m^ > 5.6 keV), we argue that the scenario in which all the DM is pro- 
duced via the DW mechanism is ruled out. We discuss however other production mechanisms 
and note that the sterile neutrino remains a viable candidate for Dark Matter, either warm or 
cold. 



1 INTRODUCTION 



A vast body of evidence points to the existence of Dark Matter 
(DM) in addition to the ordinary visible matter in the Universe. The 
evidence includes: velocity curves of galaxies in clusters and stars 
in galaxies; observations of galaxy clusters in X-rays; gravitational 
lensing data; cosmic microwave background anisotropics. While 
the DM contributes some 22% to the total energy density in the 
Universe, its properties remain largely unknown. 

The Standard Model of particle physics (SM) does not provide 
a DM candidate. The DM cannot be made out of baryons, as such 
an amount of baryonic matter cannot be generated in the frame- 
work o f an otherw ise successful scenario of Big Bang nucleosyn- 
thesis ( 15^1199 5!). In addition, current microlensing experiments 
exclude the possibility that MACHOs (massive compact halo ob- 
jects) constitu te the dominant amount of the total mass density in 
the lo cal halo dOates et al.lll995l;lLasserre et al.ll200(]|;lAlcock et al.l 
120001) . The only possible non-baryonic DM candidate in the SM 
could be the neutrino, however this possibility is ruled out by the 
present data on the large scale structure (LSS) of the Universe. 

What properties of the DM particles can be deduced from ex- 
isting observations? Some information comes from studies of struc- 
ture formation. Namely, the velocity distribution of the DM par- 
ticles at the time of structure formation affects greatly the power 
spectrum of densit y perturbations, as me asured by a variety of ex- 
periments (see e.g. iTegmark et ai]|2004 . One of the parameters, 
characterizing the influence of the DM velocity dispersion on the 
power spectrum, is the free-streaming length Afs - the distance 
traveled by the DM particle from the time when it became non- 
relativistic until today. Roughly speaking, the free-streaming length 
determines the minimal scale at which the Jeans instability can 
develop, and therefore non-trivial free-streaming implies modifi- 



cation of the spectrum of density perturbations at wave numbers 
k>A,\. 

If the DM particles have negligible velocity dispersion, they 
constitute the so-called cold DM (CDM), which forms structure in 
a "bottom-up" fashion (i.e. smaller scale objects formed first and 
then merged into the larger ones, see e.g. lPeeble s 1980). The neu- 
trino DM represents the opposite case - hot DM (HDM). Iii HDM 
scenarios, structure forms in a top-down fashion l lZerdovichll970l) . 
and the first structures to collapse have size comparable to the Hub- 
ble scale (Bisnovatyi-Kogan 1980; Bond et al. 1980; Doroshkevich 
et al. ll98ll ; iBond & Szalavlll983i) . In this scenario the galaxies do 
not have enoug h time to form, contradicting to t he existing obser- 
vations (see e.g. lWhite et al ]|l983l;|Peeble'3ll984l) . 

Warm DM (WDM) represents an intermediate case, cutting 
structure formation at some scale, with the details being dependent 
on a particular WDM model. 

Both WDM and CDM fit the LSS data equally well. The dif- 
ferences appear when one starts to analyze the details of structure 
formation for galaxy-size objects (modifications of the power spec- 
trum at momenta k > 0.5/jMpc"'). It is usually said that WDM pre- 
dicts "less power at smaller scales", meaning in particular that one 
expects smaller number of dwarf satellite galaxies and shallower 
dens ity profiles than those predicted by CDM models ( Navarro 
et al. l 19971; iKlvpin et al.lll999t lOhigna et al.ll2000l) . Thus WDM 
models can provide the way to solve the "missing satellite" prob- 
lem a nd the problem of central density peaks in galaxy-sized DM 
halos (Klvpin e t al]|l999l; iMoore et alj|l999l ; iBode et al.ll200lt 
Avila-Reese et al j200ll) . 

There exist a number of direct astrophysical observations 
which seem to contradict the N-body simulations of galaxy for- 
mations, performed in the framework of the CDM models (e.g. 
iDiemand et alj|2007l ; Istrigari et ai]|2007l) . Namely, direct mea- 
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surements of the DM density profiles in dwarf sp heroidal (dSph) 
satellites of the Milky Way favour c ored profiles jGilmore et alj 
I2OO6L l2007l : IwdlOOTt lGiimorel2007l) . ' The number of dwarf satel- 
lite galaxies, as currently observed, is still more than an order of 
magnitude below the CDM predictions, in s pite of the drastically 
improved sensitivity towards the search (see iGilmore et aT]|2007i: 
IKoposov et al.l2007 ) and resolution of numerical simulations ( Stri 



gari et al. '2007"). There seems to exist a smallest scale (~120 pc) 
at which the DM is observed (Gilmore et al. 2007: Gilmore 2007) - 
However, as of now there is no definitive statemen t about the "CDM 
substructure crisis" (see I Simon & Geha' Hmf) in regard to the 
smallest observed DM scale and [Penarrubia et al. (2007) for an 
alternative solution of the "missing satellite problem"). 

The power-spectrum of density perturbations at scales of inter- 
est for the WDM vs. CDM issue can also be studied, analyzing the 
Lyman-a forest data (absorption feature by the neutral hydrogen at 
X = 1216A at difi^ere nt red-shif ts in the distant q uasar spectra, Hui 
et al l 1997 *). This involves comparison of the observed spectra of 
Ly-or absorption lines with those obtained as a result of numerical 
simulations in various DM models. In this way one arrives at an 
upper limit on the free- streaming length of the DM particles. 

Various particle physics models provide WDM candidates. 
Possible examples include gravitinos and axinos in va rious super- 
svn mietric models (see e.g. iBaltz & Muravamij|2003l: Cembranos 
et al. I2OO6I : ISeto & Yamaguchill2007[) . Another WDM candidate 
is the sterile neutrin o with a m ass in the keV range (Dodelson & 
Widrow'[L994). Recently, this candidate received a lot of atten- 
tion. Namely, an extension of the minimal SM (MSM) with the 
three righ t-handed neutrinos was suggested ( Asaka & Shaposh- 
nikovliobs; Asaka et alj|2005t) . This extension (called vMSM) 
explains several observed phenomena beyond the MSM under the 
minimal number of assumptions. Namely, apart from the absence 
of the DM candidate, the MSM fails to explain observed neu- 
trino oscillations - the transition between neutrinos of different 
flav ors (for a review see e.g. jpogli et a Lll2006l ■ Strumia & Vis- 
sani 20061; lGiuntill2007t) . The explanation of this phenomenon is 
the existence of neutrino mass. The most natural way to provide 
this mass is to add right-handed neutrinos. Indeed, in the MSM, 
neutrinos are left-handed (all other fermions have both left-handed 
and right-handed counterparts) and strictly massless. The structure 
of the MSM dictates that right-handed neutrinos, if added to the 
theory, would not be charged with respect to any Standard Model 
interactions and interact with other matter only via mixing with 
the usual (left-handed) neutrinos (that is why right-handed neutri- 
nos are often called sterile neutrinos to distinguish them from the 
left-handed a ctive ones). Moreover, as demonstrated by Asaka & 
Shaposhnikov I2005L the parameters of added right-handed neutri- 
nos can be chosen in such a way that such a model resolves an- 
other problem of the MSM - it explains the excess of baryons over 
antibaryons in the Universe (the baryon asymmetry), while at the 
same time it does not spoil the predictions of Big Bang nucleosyn- 
thesis. For this to be true, the masses of two of these sterile neu- 
trinos should be chosen in the range 300 MeV < M2_3 < 20 GeV, 
while the mass of the third (lighter) sterile neutrino is arbitrary (as 
long as it is below Mi^j,). In particular, its mass can be in the keV 
range, providing the WDM candidate. Such a sterile neutrino can 



' For certain dSph cusped profiles are still admissible, but disfavored. Ad' 
ditional considerations rule out the p ossibility of existence of cusped prO' 



be produced in the Early Universe in the correct amount via vari- 
ous mechanisms: via non-resonant oscillations with ac tive neutri- 
nos (Dodelson & Widrow 1994: Dolgov & Hansen 2002: Abazajian 
et al. 2001; Asaka et al. 2006, 2007), via interaction with the infla- 
ton jshaposhnikov & Tkachevll2006l ). via resonant oscillations in 



the presence of lepton asymmetries Ishi & FulleJI 19991) , and have 
cosmologically long life-time. 

Finally, the sterile neutrino with mass in the keVrange would 
have other interes ting astrophy sical ap plications (see e.g. S ommer- 
Larsen & Dolgov ^2001 ) :lKusenkd fe 006); Biermann & Kusenkol 
( l200d );'Hidaka & Fulleij(l2006l);lHidaka & FuUen (i2007i): Stasielak 
et al. (i2007) and references therein). 



Existing bounds on sterile neutrino DM 

The mass of the sterile neutrino DM shou ld satisfy t he universal 
Tremain e-Gunn lower bound (Tremaine & _G unnll l979: Dalcanton 
& Hogan l200ll) : m, > 300 - 500 eV. A stronger (although model 
dependent) lower bound comes from the Lyman-o- forest analy- 
sis. Assuming a particular velocity distribution of the sterile neu- 
trino- one can obtain a relation between the DM mass and Aps and 
therefore convert an upper bound on the free-streaming length to 
a lower bound on the mass of the sterile neutrino. In the recent 
works of lSeliak et Zl j2006l) : IViel et al.l ( l2006h this bound was 
found to be 14 keV (correspondingly 10 ke V) at 95% CL in the 
Dodelson- Widrow (DW) production model dDodelson & Widrowl 
[T994) . New results from QSO tensing give similar restri ctions for 
the DW model: m, > 10 keV dMiranda & Maccidl2007l) . For dif- 
ferent models of production, the relation between the DM mass 
and the free-streaming length is different and the Lyman-a mass 
bou nd for sterile neut rinos can be as low as M, > 2.5 keV (see 
e.g. iRuchay skiy||2007l) .^ 

The sterile neutrino DM is not completely stable. In particular, 
it has a radiative decay channel into an active neutrino and a photon, 
emitting a monoenergetic photon with energy Ey = m J2 (where m, 
is the mass of the sterile neutrino). As a result, the (indirect) search 
for the DM decay line in the X-ray spectra of objects with large 
DM overdensity becomes an important way to restrict the parame- 
ters (mass and decay width) of sterile neutrino DM. During the last 



two years a number of papers appeared devoted to this task: Bo- 
yars ky et al l2006bllclldllellal; Riemer-S0rensen et al.ll2006l; Watson 



et al. I2OO6I ; iBovarskv et al.ll2007l ; lAbazaiian et al.ll2007l The cur- 



files for the Ursa Minor and Fornax jKlevna et ai]|2003albl;lGoerdt et alj 
laOOetlSanchez-Salcedo et alj2006l) . 



rent st atus of these observations is summarized, e.g., in Ruchavskivl 
( I2OO7I) . The results of the computation of st erile neutrino produc- 
tion in the early Universe dAsaka et alj2007i) . combined with these 
X-ray bounds , puts an upper bou nd on the sterile neutrino mass of 
OTj < 6 keV dAsaka et alj|2007h . This is below the lower bound 
o n the sterile neutrino DM mass from th e Lyman-o- forest analysis 
of ISeliak et al] d2006l) ; IViel et alj d2006h . Thus it would seem that 
the scenario, in which all the sterile neutrino DM is produced via 
the DW mechanism, is ruled out (the recent work bv lPalazzo et al.l 
d2007t) also explored the possibility that the sterile neutrino, pro- 
duced through DW scenario, constitutes but a fraction of DM and 

^ Sterile neutiinos are not in thermal equilibrium in the early Universe and 
therefore their velocity distribution is non-universal and depends on the 
model of production. 

^ Strictly speaking, in case of other models of production the power spec- 
trum of density fluctuations is not only characteri zed by the fr e e-stre aming 
length. Therefore, the re scaling of the results of [Seliak et alj d2006i); Viel 
et al. (2006) can be used only as the estimates and the reanalysis of the 
Lyman-dv data for the case of each model is required. 
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found this fraction t o be below 70%). However, the results of Sel- 
jak et al. (2006); Vi el et al.l (|2006|) are based on t he low-resolution 
SDSS Ly man-o; data set o ?Mc DonaId etd] l l200 ^. It was shown re- 
cen tly by Viel et al. t2007) that using high-resolution HIRES spec- 
tra teecker et al.l200'7l) one arrives at the lower limit tMj > 5.6 keV. 
Thus, the small window of masses 5.6 keV < < 6 keV remains 
open in the DW model. Therefore further improvement of X-ray 
bounds is crucial for exploring (and possibly closing) this region of 
parameters. 

It was shown in lBovarskv et al.l ( 12006 j) that the objects in the 
Local Halo (e.g. dwarf spheroidal galaxies) are the best objects in 
terms of the signal to noise ratio. The Andromeda galaxy (M31) 
is one of the nearest galaxies, excluding dwarves, that enables one 
to resolve most of its bright point sources and extract the spectrum 
of its dilfuse emi ssion. It also has a massive and well-studied dark 
matter hal o (e.g. 'Klypin et al."2002'; Wi drow & Dubinskil 120051 
iGeehan e t al. 2006; Tempel et al. 2007). The first step in such an 
analysis was done by Watson et al.l ( 12006 *) (hereafter denoted by 
W06), who analyzed the diffuse emission from the 5 central ar- 
cmin, using the data processed bv lShirev et alj j200ll) . We repeat 
the analysis of the central part of the MSI, processing more obser- 
vations, and extend the analysis to the off-centre region (5' - 13')- 
We also analyze the uncertainties in the DM distribution in the cen- 
tral part of M3 1 . The outer region of M3 1 has much fainte r diffuse 
emission than its central part (c.f. e.g. lTakahashietal.l2004 Fig. 8), 
and uncertainties in the determining of the distribution of DM in 
this region are lower. All this allows us to strengthen the restric- 
tions on the parameters of sterile neutrino DM, while using more 
conservative estimates of the DM signal. 

The paper is organized as follows. We briefly summarize the 
properties of decaying DM in Section |2] The description of DM 
in M31 and expected DM decay flux is computed in Section[3l In 
Section|4]we describe the methodology of EPIC MOS and PN data 
reduction which we perform by using two different methods: Ex- 
tended Sources Analysis Software (ESAS) and single background 
subtraction method (SBS). In Section|5]we fit the spectra and obtain 
the restrictions on sterile neutrino parameters. Finally, we discuss 
our results in Section[6l 



2 DECAYING DARK MATTER MODEL 



The flux of the DM decay 
photons s^'cm"^) is given by 

PDnir) 
4;r|DL -t- rp 



from a given direction (in 



££i f 

ins Jfm 



dr. 



(1) 



fav cone 

Here D^ is the luminosity distance between an observer and the 
centre of an observed object, poMif) is the DM density, and the 
integration is performed over the DM distribution inside the (trun- 
cated) cone - solid angle, spanned by the field of view (FoV) of the 
X-ray satellite. In case of distant objects'*, Eq. ([T) can be simplified: 



M^'"T F 
4;r£)? m/ 



(2) 



where M^^lj is the mass of DM within a telescope field of view, m, 
- mass of the sterile neutrino DM. In the case of small FoV, Eq. lO 
simplifies to 



^ Namely, if luminosity distance is mucli greater tlian tlie characteristic 
scale of the DM distribution. 



Fdm — 



where 

Sdm= I PDM(r)dr 

Jl.o.s. 



(3) 



(4) 



is the DM column density (the integral goes along the line of sight), 
£1 <K 1 - FoV solid angle. 

The decay rate of the sterile neutrino DM is equal to ( Pal & 
Wolfenstein 1982; B arger et al. II 19951) ^ 



9aGl 
I024;r'' 



sin^(26')m; =a 1.38 ■ 10"^°5" 
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1 keVJ 



(5) 



Here /n, is the sterile neutrino mass, 6 - mixing angle between ster- 
ile and active neutrinos. From a compact cloud of sterile neutrino 
DM we therefore obtain the flux: 



6.38 • 10^ 



keV 



kpc 



D, 



sin- (26) 



I keVJ 



.(6) 



3 ANDROMEDA GALAXY (M31) 

M31, or Andromeda galaxy, is one of the nearest galaxies, exclud- 
ing dwarves; it is located a t the distance = 784 ± 13 ± 17 kpc 
dStanek & Gamavicljl99"8l) . Its proximity allows us to resolve most 
of its point sources and extract the spectrum of diffuse emission of 
its central part. 

Available XMM-Newton jjansen et al.ll200ll) observations 
cover the region of central 15' of M31 with exposure time greater 
than 100 ksec (see TableO. W06 used the XMM-JVewton data on 
centra l 5' of M3 1 (observation Q 11257(8401 pr ocessed by Shirey 
et al. (2001'), exposure time about 30 ksec) to produce restrictions 
on the parameters of sterile neutrino DM. The sufficient increase 
of photon statistics enables us to analyze the outer (5' -13') faint 
part of M3I, which, however, has a significant mass of DM (see 
Section l3.1l below). 

In this work we will analyze two different spatial regions of 
Andromeda galaxy: region circles, which corresponds to 5' cir- 
cle around the centre of M31, and region ring5-13, which corre- 
sponds to the ring with inner and outer radii of 5' and 13', respec- 
tively. 



3.1 Calculation of DM mass 

To obtain the restriction on parameters of the decaying DM, we 
should calculate the total DM mass M^^, which corresponds to 
both spatial regions: circles and ring5-13, both with and with- 
out resolved point sources. To estimate the systematic uncertainties 
of the evaluation of the DM decay signal and to find the most con- 
servative estimate for it, we analyze va rious avai lable DM profiles 
avpin e t al. 200 21; Iwidrow^ & Dubinski 200j; 



i Kerins etalJ 2001 



iGeehan et alhooe 



Carignan et all 20061; iTempel et aU2007l) : 



^ Our decay rate is 2 times smaller than the one used in W06. This is due to 
the MaioraiMjiature_ofthe_sterileji^^ 

et al. 1995). The final constraints for a Dirac particle would thus be 2 times 
stronger. 
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Obs. ID 


Starting time, UTC 


Filter 


Cleaned MOS1/MOS2/PN exposure, ks 


0112570401 


2000-06-25 08:12:41 


Medium 


30.8/31.0/27.6 


0109270101 


2001-06-29 06:15:17 


Medium 


40.1/41.9/47.4 


0112570101 


2002-01-06 18:00:56 


Thin 


63.0/63.0/55.3 



Table 1. Observations of the central part of M3 1 , used in our analysis. 



Model 


circle5 


ring5-13 


13 arcmin sphere, MC result 


13 arcmin sphere, analytical result 


Kl, with sources 


3.27 ±0.01 


12.49 ± 0.03 


5.84 ± 0.02 


5.84 


K2, with sources 


11.88 ±0.03 


23.75 ± 0.09 


20.76 ± 0.09 




GFBG, with sources 


6.59 ± 0.02 


20.46 ± 0.06 


13.40 ± 0.03 


13.39 


KING, with sources 


6.68 ± 0.01 


24.61 ±0.05 


14.80 ± 0.02 


14.80 


MOORE, with sources 


7.34 ± 0.02 


19.48 ± 0.02 


13.79 ± 0.02 


13.78 


N04, with sources 


7.68 ± 0.03 


22.89 ± 0.07 


15.16 ±0.06 


15.18 


NFW, with sources 


11.08 ±0.04 


40.5 ±0.1 


22.3 ±0.1 


22.25 


BURK, with sources 


6.71 ± 0.02 


27.97 ± 0.03 


15.90 ±0.05 


15.90 


KER, with sources 


5.35 ± 0.02 


22.45 ± 0.04 


11.56 ±0.03 


11.56 


M31A, with sources 


5.95 ±0.01 


16.45 ± 0.02 


11.03 ±0.02 




M31B, with sources 


4.99 ± 0.01 


14.24 ±0.01 


9.40 ± 0.02 




M31C, with sources 


5.60 ±0.01 


16.12 ±0.01 


10.29 ± 0.02 





Table 2. DM mass (in lO' M©) inside regions, used in ou r analysis: res ults of our Monte Carlo integration. The point s ources are not exclud ed h ere. The 95% 
statistical errors are also shown. The DM distributions of iKlvpin et al.l ( 2002) (before and after ad iabatic contraction), Geehan et al.l i200d) and lKerins et alj 
)200l[) are marked as "Kl", "K2", "GFBG" and "KER", respectiv ely. The DM dis tributions from lTempel et alj )2007l) are marked as "KING", "MOORE", 
"N04", "NFW" and "BURK" (see text). The DM distributions from lwidrow & Dubinski (2002) are marked as "M31A", "M31B" and "M31C". 



• (Kl) Before^ adiabatic contraction stage. lKlvpin et al] ( |2002|) as- 
sume that DM distri bution is purely Navarro-Frenk- White (NFW) 
jNavarro et aLlll997l) : 

PDM(r) 



An [log(l +C)- C/(l + O] r{r + 



(7) 



The paramet ers of this NFW di stribution (in terms of the favored 
CI model of iKlvpin et ai]|2002h are: M,,> = 1.60 x 10'^ Mq\ r, = 
25.0 kpc; C = 12. 

• (K2) This non-analytical model is the result of adiabatic contrac- 
tion of the Kl profile, described above. To obtain it, we extract the 
data from the Fig. 4 of iKlvpin et al.l i2002l) . In the top part of this 
figure the dot-dashed curve is the contributio n of the DM halo t o 
the total M31 mass distribution (CI model of ' Klvpin et al.ll2002l) . 
As the precise form of this mass distribution is not analytic, we 
scanned this curve and produced the file with numerical values of 
enclosed mass Mouir) within the sphere of radius r. After that, we 
interpolated the Mouir), and evaluated the radial density distribu- 
tion 

1 dMpuir) 

• (GFBG) Preferre d Navarro-Frenk- White distribution from Gee- 
han et al. ( 2006.) : M,„> = 6.80 x 10" Ms; = 8.18 kpc ; C = 22. 

• (KER) Isothermal profile used in lKerins et alj ( l200lh : 



PKER(r) ■ 



r > R,„ 



(9) 



where pi,(0) = 0.23Mo pc = 2 kpc, R^^^. = 200 kp c. 
• (M31A-C) Profiles of i Widrow & Dubinskil ( |2005|) . In this pa- 
per the authors propose several models, which differ by the relative 
disk/halo contribution. These non-analytical models (M31a-d) in- 
corporate an exponential disk, a Hemquist model bulge, an NFW 



halo (before contraction) and a central supermassive black hole. 
The stability against the formation of bars was numerically stud- 
ied.' 

We also use density distributions from the recent paper of Tempel 
et al. (200i). The main aim of this paper is to derive the DM den- 
sity distribution in the central part of M31 (0.02-35 kpc from the 
centre). 



• (K ING) Modified isothermal profile jKindl 19621 ; lEinasto et all 
fl97i) : 



Pisoir) ■ 



r> ro. 



1 + 4 



r < To, 



(10) 



where po = 0.413Mo pc 



-3 ,. _ 



1.47 kpc, rp = 1 17 kpc. 



(MOORE) Moore profile jMoore et alJl 19991) : 

Pc 



Pmoore(i') ■ 



(11) 



where p, = 4 A3 ■ IO^'Mq pc ^ = 17.9 kpc 

• (N04) Density distribution of lNavarro et alj ( |2004|) : 

PNQ4(r) = Pc exp 



-ip-1 



(12) 



where parameter a, accordi ng to simulations, equals to 0.172 
0.032 jNavarro et ai]|2004l) . For N04 we take a = 0.17, p, 
6.42 ■ lQ-^Mapc-\ r, = 11.6 kpc. 



(NFW) Navarro-Frenk- White profile: 

Pc 



PNEwif) ■ 



(13) 



* In contrast to the other models, this model does not describe the cunent 
DM distribution, but helps our understanding the time evolution of DM 
mass inside constant FoV. 



' We do no^^seJhe^ourtl^Tiode^_^M31d]j_^ecausejn Widrow & Dubin- 
ski (2005) it was found that this model develops a bar, which rules it out 
experimentally. 
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Model 


circles 


Removed from circles , % 


ring5-13 


Removed from ring5-13, % 


Kl, without sources 


U. /O / ± U.UU4 


/O.D 


y. / 1 ± U.Uz 


22.3 


K2, without sources 


oil 1 r\ Ai 

2.31 ± U.U2 


80.4 


1 o nci 1 rv AO 

18.09 ± 0.08 


23.9 


GFBG, without sources 


1.48 ± U.Ui 


/ /.4 


i J. / / ± U.UO 


Z.J.U 


KING, without sources 


1.64 ±0.01 


75.5 


18.99 ±0.06 


22.9 


MOORE, without sources 


1.52 ±0.01 


79.2 


14.98 ± 0.03 


23.1 


N04, without sources 


1.70 ±0.02 


77.7 


17.62 ± 0.05 


23.0 


NFW, without sources 


2.59 ± 0.01 


76.7 


31.34 ±0.07 


22.5 


BURK, without sources 


1.67 ±0.02 


75.1 


21.68 ±0.02 


22.5 


KER, without sources 


1.33 ±0.01 


75.0 


17.42 ± 0.04 


22.5 


M31A, without sources 


1.24 ±0.01 


79.3 


12.66 ± 0.02 


22.9 


M31B, without sources 


1.04 ±0.01 


79.1 


10.98 ±0.01 


23.0 


M31C, without sources 


1.21 ±0.01 


78.4 


12.43 ± 0.01 


22.9 



Table 3. DM mass (in lO' Mq) without point sources: results of our Monte Carlo integration. The fraction of DM, removed together with the point sources, 
is also shown. All notations are the same as in previous table. 




Figure 1. Selected regions in the central part of M3 1 (shown in linear scale). 
Small circles correspond to excluded point source regions, large circles have 
radius of 5 and 1 3 arcmin. 
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Figure 2. M3 1 DM column density versus off-centre angle as result of our 
Monte Carlo integration, based on DM profiles of Sec. 13. II (Point sources 
are not excluded). 



where = 5.20 ■ 10"- Mq p c~\ = 8.31 kpc. 
• (BURK) Burkert profile teurkertll 19951) : 



Pburk{^) 



Po 



(14) 



(l-t)(l-f)' 
where po = 0.335Mo pc"^, = 3.43 kpc. 

The computed DM masses within the FoV for all these pro- 
files are shown in Table [2] We see that for the model used by 
W06 (model K2 in our notations), our estimate of the DM mass 
within the central 5' coincides with the value used in W06: M5 = 
(1.3±0.2)- IO'^Mq. Notice, however, that to obtain the diffuse spec- 
trum, we extracted all point sources, resolved with the significance 
> 4<T. Each source was removed with the circle of the radius of 
36" (see Sec. l4. ll for details). This led to the reduction of the area 
of the FoV by about 70% in case of circles region (c.f. Fig.[T](. 
As the density of the DM changes with the off-centre distance and 
this change can be significant (c.f. Fig.|2j, we performed the inte- 
gration of the DM density distribution over the FoV with excluded 
point sources. To calculate the DM mass in such "swiss cheese" 
regions (Fig.lTJ, we used Monte Carlo integration. The results are 
summarized in the Table [3] 

To check possible systematic effects of our Monte Carlo in- 
tegration method, we also obtained the values of enclosed mass 



inside the 13 arcmin sphere, and compared them with analytical 
calculations (wherever possible). Such an error does not exceed the 
purely statistical error of numerical integration (see Table|2ll. 

As one can see from Tables [2l3l the most conservative DM 
model, d escribing regions circles and ring S- 13, is the model 
M31B of IWidrow & Dubinskil ( |2005|) . Therefore, to obtain restric- 
tions on the DM parameters in what follows, we will use the DM 
mass estimates based on this model. 

For the DM distributions listed above, we also build the DM 
column density 5 dm (given by Eq. l|4j) versus off-centre angle. The 
result is shown on Fig. [2] It is clearly seen that, in the off-centre 
regions, there is still a lot of DM, and, together with the fact that the 
surface brightnes s of X-ray diff^use em ission falls rapidly outside 
the central 5' (c.f. lTakahashi et alj2004 ). improving the restrictions 
of W06 by analyzing the off-centre 5'-13' ring. Moreover, as one 
can see from Table [3] and Fig. |2] the uncertainty of DM in this 
region is less than in the circles region. 

To estimate the additional contribution from the Milky Way 
DM halo in the di rection of M31, we use an is othermal DM dis- 
tribution (as e.g. in iBovarskv et al.ll2006dll2007h . The DM column 
density is equal to 



■K((t>), 



(15) 
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where 
model, Tq 
and 



170 km s- 



4 kpc - parameters of isothermal 



8 kpc - distance from Earth to the Galactic Centre, 
f +arctan[i^], cos^>0 



(16) 



arctan 



r© cos J 



cos < 0. 



Here is defined via cos 4> = cos / cos b for an object with galactic 
+ sin" 0) . For Andromeda galaxy 
(/= 121.17°, i = -21.57°, i.e. 0= 118.77°) one obtains 



Smw.dm ~ 6.2 • 10 • cm ^ 



3.5 X 10^' keV • cm"- 



(17) 



According to Fig. ID the MW contributes < 5% to the total DM 
column density along the central part of Andromeda galaxy, and 
therefore will be neglected in what follows. 
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4 DATA REDUCTION AND BACKGROUND 
SUBTRACTION 

To obtain restrictions on the parameters of the sterile neutrino, we 
need to analyze diffuse emission from faint extended regions of 
M31. There exist several well-developed background subtraction 
procedures for the difi^use sources (see, for instance, XMM-JVewto n 
SAS User Guide MNevalainen et al.ll200llRead & PonmarJl2003h . 
In this paper we use two methods of background subtraction: 

4.1 Extended Sources Analysis Software (ESAS) 

This method, recently developed by ESAC/GSFC team', allows 
one to subtract instrumental and cosmic backgrounds separately. It 
seems to be better than the subtraction of the scaled blank-sky back- 
ground, averaged through the entire XMM-Newton Field of View 
(see next subsection for details), as instrumental and cosmic back- 
grounds (due to their different origin) have different vignetting cor- 
rection factors. ESAS models instrumental background from "first 
principles", using filter-wheel closed data and data from the unex- 
posed corners of archived observations. Using this software, we are 
assured that no DM line can be in our background, in contrast with 
the "black sky" background subtraction method and, espe cially, lo- 
cal background subtraction (used e.g. in IShirevetal.l ('2Q01) to pro- 
duce the diffuse spectrum of central 5' of M31). The price to pay is 
the necessity of modelling cosm ic background. 

To prepare the EPIC MOS jTurner et a0200ll) event lists, we 
used the ESAS script mos-filter. After running mos-filter, 
we produced cleaned MOS images in sky coordinates, which were 
used to obtain the mosaic image (with the help of SAS v.7.0.0 tool 
emosaic). We used these event lists and images to find the point 
sources using SAS task edetect.chain. Source detections were 
accepted with likelihood values above 10 (about 4cr). We found 
243 point sources in this way. After that, we excluded each of them 
within the circular region of the radius 36", which corresponds to 
the removal of ~ 70 - 85% of total encircled energy, depending 
on the on-axis angle (see XMM users handbook'" for details). The 



http: //xmm. esac.esa. int/external/xmm_user_ 
support/documentation/sas_usg/USG 
' We use ESAS version 1.0. 

'" http : //xHim . esac . esa . int/external/xrmti_user_ 
support/documentation/uhb 



Figure 3. Observed spectrum (top) and modelled instrumental background 
(bottom) MOSl from ObsID 01 12570101, region ring5-13. It can be seen 
that the spectrum and modelled background almost coincide for E > 7keV. 

constructed mosaic image with detected point sources and selected 
regions is shown in Fig.[T] 

We obtained the MOSl and MOS2 spectra and constructed 
the corresponding background with the help of ESAS scripts 
mos-spectra" and xmm-back, respectively. 

Finally, we grouped the spectra with corresponding response 
and background files with the help of FTOOL grppha, a part of 
HEASOFT v6.1. To ensure Gaussian statistics, the minimum num- 
ber of counts per bin was set to be 50. 

The ESAS method of background subtraction, however, has 
several difficulties. The number of fitting parameters substantially 
increases, hence it is harder to find true minimum of The quan- 
titative analysis of the 1.3-1.8 keV energy range is also not possi- 
ble, because of the presence of t wo strong unmodell ed instrumental 
lines (see Figs.[3l|3l. EPIC-PN (Striider et alj|200lb data reduction 
is not yet implemented in ESAS. Therefore, to cross-check the re- 
sults obtained with the help of ESAS software, we also processed 
EPIC d ata with the help of the blank-sky data subtraction (SBS) 
method dRead & PonmaiJl2003h . 



4.2 Blank-sky background subtraction (SBS) 

We processed the same M31 observations (Table[T) as in the previ- 
ous Section, using both MOS and PN data. To subtract the blank- 
sky background we firstly cast it at the position of M3 1 with the 
help of the script skycast '-, written by the XMM-Newton group 
in Birmingham. The scaling coefficient was derived by comparing 
count rates for £ > 10 keV from source regions and background 
sample. To produce spectra, ARE, RMF and to group them correctly 
(we needed to extract them from non-circular regions), we modi- 
fied the Birmingham script createspectra.''' The spatial regions 
were chosen similarly to those in Sec. 14. II so it would be possible 
to compare the results of the two different methods (see Sec l5.3t . 



To produce correct RMF file, we changed in the script mos-spectra 
option rmfgen detmaptype=psf to rmfgen detmaptype=dataset. 
http : //www . sr . bham . ac . uk/xirmi3/skycast 
ht tp : //www ■ sr . bham . ac . uk/ximii3/createspectra| 
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Figure 4. Folded spectra from ring 5 - 1 3 region (by ES AS method), with 
excluded point sources. The presence of two unsubtracted instrumental lines 
at 1.49 keVand 1.75 keVis clearly shown. 
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Figure 5. Unfolded spectra and best-fit model from ring5-13 region (by 
ESAS method), with excluded point sources. The "line forest" at energies 
lower 2.0 keV is clearly visible. 

When analyzing PN data, we found that the role of out-of- 
time (OOT) events was significant. This is due to the fact that 
the rate of OOT events is proportional to the total rate inside the 
full PN FoV rather than the rate of diffuse emission (outside ex- 
cluded point sources). Therefore, it was necessary to remove the 
OOT events from the PN event lists. Most of the OOT events (from 
the bright point sources) form strips in the images and can be easily 
removed with the help of spatial filtering. This additional filtering 
also slightly reduced the possible DM signal, which was (in this 
outer region) nearly proportional to BACKSCALE keyword. This was 
accounted for when producing SBS PN restrictions. 



Notice that the fit results obtained by two background subtraction 
methods (ESAS and SBS) coincide within the 90% confidence in- 
terval (Table |4l.''' Also shown in Table |4l are t he results of Taka- 
hashi et al. (2004), who analyzed diffuse emission in the central 
6' of M31.'^ Below we discuss separately the fitting of ESAS and 
SBS spectra. 



5.1 ESAS spectra 

We build 0.5 - 10.0 keV MOS spectra of circles and ring5-13 
regions for 3 observations from Table [T]'* Thus for each spatial 
region we have 6 spectra to fit - from observations with MOS 1 and 
M0S2 cameras. We fix the model parameters to be equal for all six 
spectra from the same spatial region (except for normalization of 
the remaining soft proton background, as the spectra from different 
observations are slightly different). 

Since ESAS software subtracts only the instrumental back- 
ground component, the remaining cosmic background should be 
modelled. The cosmic background component is modelled with 
the help of Xspec model apec+(apec+pow) '-wabs, according to 
the ESA S manual. A cool (~ 0.1 keV). unabsorbed apec ( Smith 
et al. | 200l[ ) component represents the thermal emission from the 
Local Hot Bubble. The hot (~ 0.25 keV), absorbed apec compo- 
nent represents emission from the hotter halo and/or intergalactic 
medium. The last, absorbed powerlaw component with powerlaw 
index T = 1.41 represents the unresolved background from cos- 
mological sources. We kept its normalization fixed for each region; 
it corresponds to 8.88 ■ 10"^ Xspec units per square arcmin, or to 
10.5 photons keV"' s"' cm"^ sr"'. The corresponding hydrogen 
column density in w abs was left to vary below its Galactic value 
iiH = 6.7 • iO^°cm-^ tMorrison & McCaimnorilll983l) . To model the 
soft proton contamination, we used bknpow/b model (we fix its 
break energy at 3.3 keV), where index /b means that this compo- 
nent is not folded through the instrumental effective area (in Xspec 
versions 11 and earlier). 

The diskbb+bb odv (the same as the LMXB model in Taka- 
hashi et al. 2004) component describes the point sources, which 
were not excluded. We fitted the dilTuse M31 c omponent in outer 
regions with the help of the sum of three vmekal dMewe et al J 1 986t 
iLiedahl et aLlll995h models with fixed temperatures and abun- 
dances. The wabs column density was fixed at its Galactic value. 



5.2 SBS spectra 

We fitted the data from MOS and PN cameras, pro- 
cessed using SBS method (both separately and combined). 
As both cosmic and instrumental background is sub- 
tracted in SBS method, we fitted MOS and PN spectra on 
wabs"- (diskbb+bbody+vmekal+vmekal+vmekal) Xspec model 
at the energy range 0.6-10.0 keV (0.6-12.0 keV in case of PN 
camera). The reduced obtained by fitting our spectra are shown 
in Table|5] fit parameters are shown in Table|4] 



5 FITTING THE SPECTRA IN XSPEC AND PRODUCING 
RESTRICTIONS 

After we have prepared the data (with ESAS and SBS back- 
ground subtraction methods) we fitted obtained spectra with real- 
istic model (u sing Xspec spectral fitting package version 11.3.2, 
lAmaud|[T993) . The results of our fits are shown in Tables |4l |5] [6] 



The value of normj^ also coincides within 90% confidence interval if 
one propagates the uncertainty of blank-sky background normalization. 

The appreciabledifference between ouremrsandthoseof Takahashi 
et al. (2004) is due to the fact that we did not fix the metal abundances 
equal to each other. This was essential for our purposes, because of the 
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Table 4. Model parameters from regions circles and ring5-13. Also shown are 90% confidence intervals for fitted parameters. Results of lTakahashi et alj 
(6' circular region in this case) are marked as "TOKM". 
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Figure 6. Folded MOSl spectra from circles region, ObslD 01 12570401, 
with (top) and without (bottom) point sources. 
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Figure 7. Folded spectra and best-fit model from circles region, with 
excluded point sources. 



Region 


Reduced x'' 


Number d.o.f. 


ESAS, circles 


1.071 


399 


SBSMOS, circles 


1.102 


371 


ESAS, rii\g5-13 


1.109 


1608 


SBS MOS, ringS-13 


0.994 


1735 


SBS PN, ringS-13 


1.007 


2754 


SBS PN-OOT, rings- 13 


0.995 


2715 


SBS MOSPN-OOT, ringS-13 


1.009 


4082 



Table 5. Reduced for our regions. 



5.3 Producing restrictions on sterile neutrino parameters 

In this subsection we describe two different techniques of search- 
ing for the narrow (compared to the spectral resolution of XMM- 
Newton) decay line in the spectra, processed by ESAS and SBS 
methods. 

As shown on Fig. [5] above 2.0 keV there are few emission 
lines in the model of the spectrum of M31, and continuum emis- 
sion dominates. In this ca se, it is possible to apply the "statistical" 
method, discussed e.g. in lBovarskv et al.l ( l2006dl . Namely, after 
fitting the spectra with the selected models (Sees. l5T[j5. 21 above), 
we add an extra Gaussian line with the help of Xspec command 
addcomp. We then freeze its energy Ey, leave the line width cr to 
vary within 0-10 eV, and repeat the fit. For each line energy, we re- 
fit the model and derive an upper limit on the flux in the Gaussian 
line, allowing all other model parameters to vary. In particular we 
allow the abundances of heavy elements, that produce the thermal 
emission lines to vary. This produces the most conservative restric- 
tions as the added line could account for some of the flux from 
the thermal components. After that we calculate the 3(T error with 
the help of Xspec command error {line norm) 9.Q. To obtain 
conservative upper limits, we allow as much freedom as possible 
for the parameters of the thermal model. The 3cr upper limit on 
the DM line flux is shown in Fig. [8] These flux restrictions can be 
turned into constraints on parameters of the sterile neutrino (m , and 
s'm'{29)), using Eq. ^ and the value of the M-^^ from the Table[3] 
for the model M31B. 

Below 2.0 keV, there are a lot of strong emission lines, which 
dominate over the continuum, creating a "line forest". As the in- 
trinsic widths of these lines are much more narrow than the spectral 
resolution of EPIC cameras of XMM-Newton, and the abundances 
of various elements are known with large uncertainties, it is very 
hard to reliably distinguish these emission lines from a possible 
DM decay line. Therefore, to produce robust constraints, we apply 
the "full flux" method below 2 keV. In this method, we equate the 
DM line flux to the full flux plus 3 flux uncertainties over the energy 
interval AE equal to the spectral resolution of the instrument." 

We also produce model-dependent "statistical" constraints be- 
low 2.0 keV. To reduce model uncertainty, we fix most metal abun- 
dances at their va lues known from optical observations of M31 (Ja- 
coby & Ciardullo Tl999l : |jacobv & Fordll 19861; iDennefeld & Kuntlj 



clear presence of the "line forest" at energies below 2.0 keV (see Sec. |5.3| 
and Fig.|5). 

" We exclude the region 1.3-1.8 keV due to the presence of two strong 
unmodelled instrumental lines, see Sec. 14. II 

" To find the proper value of &E, we fold thin Gaussian line with ap- 
propriate RMF, and then evaluate FWHM of obtained broadened line. The 
FWHM AE. calculated in such a manner, slowly increases with line energy 
and changes from 0.18 keV to 0.21 keV in the 0.5 - 2.0 keV energy region. 
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Figure 8. 3(T upper limit on the DM line flux (the region of parameter space above the curves is excluded). Left panel: upper Hmits from the different spatial 
regions for the spectra, processed by ESAS method. Right panel: upper limits for the ringS - 13 region for both ESAS and SBS methods. 
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0.2..1.7 
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0.3..2.7 


0.2..0.5 



Table 6. Abundances from optical observations (in solar units). Our allowed range of abundances, used for construction the model-dependent restriction (see 
Sec. l5.3t . is also shown. 



ll98ll:lBlairetal.ll982h . The confidence ranges of these abundances 
are shown in Table |6] 

To compare our results with previous work on M3 1 (Watson 
et al. l200d hereafter W06) we performed full flux analysis in the 
whole region of energies of the MOS camera of XMM-Newton. 
The results are shown in Fig.[9] One can see that our full flux results 
from circles region are somewhat weaker that the corresponding 
results of W06 (by a factor 2-3 in the region m, ~ 4 keV; more 
than an order of magnitude at hj, < 2keV and > 12keV). There 
are several reasons for this. As discussed in Sec. l3.1l we use an ~ 8 
times lower estimate for the DM mass within the FoV, because we 
use the mo re^recen t and more conservative DM profile of Widrow 
& Dubinski ( 120051) and compute the amount of DM by explicit inte- 
gration over the FoV with removed point sources. At the same time, 
comparing our diffuse spectrum (Figs. [617] ) with Fig. 1 in W06, we 
see that the intensity of our diffuse spectrum is ~ 2 - 3 times lower 
(due to the ~ 4 times larger number of point sources removed). 
Therefore, one would expect a factor 2-3 difi^erence between our 
results (as indeed is seen at m., ~ 4 keV). 

An additional discrepancy at low energies is due to the differ- 
ent choice of the energy bin intervals. In W06 the energy bin inter- 
val was chosen according to the empirical formula A£ = Ey/30 = 
mJ60, while we have determined it using the XMM-Newton re- 
sponse matrices (as described in footnote [T7] above). The difference 
is most prominent at low energies: e.g. at £ ~ 1 keV we obtain 
AE » 0.2 keV, which is ~ 6 times bigger than the value, used 
by W06. Therefore, at small energies we would expect constraints 
about an order of magnitude lower than those of W06, as Fig. |9] 
indeed demonstrates. 

The other important effect, seen in Fig. [9l is the high-energy 
behaviour. Our restrictions remain nearly constant for ra, > 12keV 
(Ey > 6 keV), in contrast to the steeply decreasing results of W06. 



This is due to the fact that W06 used an energy-averaged count- 
rate-to-flux conversion factor (i.e., the telescope effective area): see 
Sec. IV of W06. However, the effective area of the XMM-Newton 
MOS cameras declines sharply with energy, essentially going to 
zero at 9-10 keV.'^ Therefore, after a proper conversion, a constant 
count rate at high energies, assumed by W06 would correspond to a 
sharply rising physical flux in photons/(s • cm'), which is of course 
incorrect. We performed a full data analysis, taking into account the 
dependence of the effective area on the energy and our constraints 
weaken sharply at high energies. This effect is well-known and is 
present in many papers that perform spectral analysis of XMM- 
Newton or Chandra data. 

Our final constraints are shown in Fig. [TO] At masses m, > 
4 keV (energies Ey > 2 keV) we use the results of statistical con- 
straints from the ringS - 13 region. To produce the final restriction, 
we choose, for each value of m„ the minimal value of sin'(2(?). For 
m, < 4keV (Ey < 2 keV) we plot both the model-independent (full 
flux) and the model-dependent constraints. The restrictions of Bo- 
yarsky et al. 1 12007,) and .Watson et aL (.2006.) are shown for com- 
parison. 

The high-energy behaviour of our final statistical constraints 
differs from that of in Fig. |9] There are several reasons for this. 
Firstly, in Fig.|9]we showed the full flux restrictions from the MOS 
camera (to compare our results with those of W06), while in Fig. 1 101 
we used the combined constraints from both MOS and PN cam- 
eras. The PN camera has a wider energy range: its effective area 
decreases only above £ » 10 keV", which explains the weak- 

18 p(jj. pj^ camera this happens at ~ 12 keV (c.f. Fig. [8). 
" XMM-Newton Users Handbook, Sec. 3.2.2.1, 

http: //xmm. esac . esa. int/external/xmmjuser. 
support/docuinentation/uhb_2 . 5 
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Figure 9. Our limits on ^m ,, sin^(20)^ parameters, obtained by using the full 
flux method from different spatial regions of M3 1 (a region of parameter 
space above a curve is excluded). The restriction from W06 is shown for 
comparison. 
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Figure 10. Restrictions on (m,, sin (26)) plane. The strongest previous lim- 
its of Bovarskv et al. ( 2007) as well as results of W06 are shown for com- 
parison. The region above the curve is excluded. 



ening of constraints on Fig[TO] for m, > 20 keV. The "peak" at 
» 16 - 18 keV, is due to the presence of strong Cu instrum en- 
tal Hnes in the PN background spectrum dStriider et alj|200ll see 
also Fig.[8]l. This region has, thus, higher errors, which weaken the 
constraints. Finally, we used several jointly fitted spectra (up to 9 
in MOSPN-OOT dataset) in our "statistical" method, as opposed 
to the restrictions in Fig|9] where we used only one spectrum. The 
combination of several spectra improves the bounds, as statistical 
errors decrease. 



out to be the most conservative among those which studied the DM 
distribution in the inner part of MSl.^" 

We found that exclusion of numerous point sources from the 
central part significantly improves our limits, therefore we have 
also calculated the DM mass in such "cheesed" regions with the 
help of Monte Carlo integration. 

As the surface brightness is low in the selected regions, the 
choice of the background subtraction method is important. We pro- 
cessed XMM-Newton data from these regions with the help of 
two different background subtraction techniques - the Extended 
Sources Analysis Software (ESAS), and the blank-sky background 
subtraction (SB S), using the blank-sky background dataset o f Read 
& Ponman (2003). We have shown that these totally different back- 
ground subtraction methods give similar results. 

To compare our results with the previous work on M3 1 ( Wat- 
son et al. 12006", W06), we obtained the full flux restriction from the 
central 5' of M31. Our full flux results (shown in Fig.|9]l mostly 
reproduce the results of W06, up to differences arising from our 
more conservative estimate of expected DM signal and proper data 
analysis (see Sec. l5.3l for detailed discussion). 

Our final upper limits (both model-dependent and model- 
independent) are shown in Fig. [TO] We improved the previous 
bounds of W06 on sin^(20) by as much as an order of magnitude 
for masses 4 keV < m, < 8 keV. Due to the significant low-energy 
thermal component in M31 diffuse emission, to produce the model- 
independent constraints, we have used the "full flux" method for 
m, < 4.0 keV (i.e. Ey < 2.0 keV). In this regio n, the strongest 
constraints remain those of iBovarskv et alj ( l2007h . We have also 
produced model-dependent constraints for Ey < 2.0 keV, using the 
"statistical" method; in this case we found the best-fit model by 
fixing the metallic abundances at the level of optical observations. 

The comparison of our upper li mit with the lower bound on 
sterile neutrino pulsar kick mechanism jFuller et alj2003l) improves 
the previous bounds and can exclude part of the parameter region 
(for 4 keV< m, < 20 keV). 

Finally, it should be noticed that although throughout this pa- 
per we were writing about the sterile neutrino DM, the results of 
this work are equally applicable to any decaying DM candidate 
(e.g. gravitino), emitting photon of energy Ey and having decay 
width F. Our final results in this case are presented in Fig. [TT] 
For other works discu ssing cosmological and astrophysical effects 
of d ecaying DM see de Ruiula & Glashowl I1980I): Berezhian i 



et al. (ll987l):lDoroshkeyich et alJ l ll989l) : lBerezhiani et al.l ( ll990h : 
iBerezhiani & Khlopovl (1 1990^. An extensive review of the results 
can also be found in the book bv .Khlopov, C1997.) . 



6.1 Sterile neutrino in Dodelson-Widrow model 

The results of this work have important consequences to one of the 
production models for the sterile neutrino, the so-called "Dodelson- 
Widrow" (DW) scenario - produ ction through (non-resonan t) os- 
cillations with an active neutrino jDodelson & Widrowlll994 . The 



6 RESULTS AND CONCLUSIONS 

Using available XMM-Newton data on the central region of the 
Andromeda galaxy (M31), we obtained new restrictions on sterile 
neutrino Dark Matter parameters. We analyzed various DM distri- 
butions for the central part of M31, and obtained a conservative 
estimate of the DM mass inside the central 13', using the model 
M31B of lWidrow & Dubinskil fcoOSh . This DM distribution turned 



We would like to notice, however, that in the work lKerinsH2004 . a num- 
ber of "extreme" (i.e. maximizing contributions of disk, spheroid or halo) 
models are considered. Some of these models would reduce an estimated 
DM signal from the inner 13' (and correspondingly our limits) by a factor 
~ 2. We chose to use the family of models, shown on Fig.|2] as they qualita- 
tively agree with each other and do not contain any "extreme" assumptions. 
However, below, in deriving a model-dependent upper limit of the mass of 
the DM particle, we will introduce an additional penalty factor, to account 
for this and other possible systematic uncertainties. 
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Figure 11. Constraints on the decay wid th F of any radiatively decaying 
DM from this work (marked "M31") and lBovarskv et al.l j2007h (marked 
"MW"). The shaded region of parameters is excluded. 




Figure 12. Current X-ray constraints, combined with the DW production 
model. Colored regions are excluded. The grey region shows the range of 
param eters which give correct abundance in the DW model iAsaka et alj 
I2OO7I) . The color_^iade^j|egionsjii^kdieresri 

yarsky et al. l2006j) . "MW" teovarskv et al.ll2007h and "M31" (this work). 
Model-dependent restrictions from M31 for m, < 2 keV are shown in 
(green) dashed line. 



computation of the abundance is complicated in this case by the 
fact that the production mainly happens around the QCD transition 
and therefore QCD contributions are hard to compute (see Asaka 
et al. l2006[ and refs. therein). A first-principles computation, taking 
in to account al l QCD contributions in a proper way, was performed 
in lAsakaetal] ( l2007l) . 

We compare the results of this computation with X-ray bounds 
obtained in this work and previous works in Fig.[T2] The upper and 
lower dashed lines, bounding the grey area, correspond to the DW 
production scenario when all hadronic uncertainties are pushed in 
one or another direction; the thick central line corresponds to the 
most probable relation between m, and sm^(2d). Upon comparison 
with X-ray bounds, we find that the upper bound on the DM mass in 
the DW scenario is reliably below m, < 4 keV (even if we push our 
X-ray bounds up by a factor of 2, to account for some yet unknown 
systematics and push all the uncertainties in hadronic contributions 
to the DW production in one direction). 



s imp: 

lAsaka et alj ( l2007h . Notice that other bo unds on /Wj, that ap peared 
in the literature (e.g. < 3.5 keV of IWatson et alj ( 1200 6) and 
Ills < 3 ke V of iBovarskv et al.l ( l2006d) ) were based on the com- 
putations of lAbazaiianl 1 20061) . which did not take into account all 
QCD contributions. 

O ur present results may be combined wit h the Lyman - or anal - 
ysis of ISeliak et all j2006l) : Iviel et al.1 {2W6); Viel et al. (2007|). 
As follows from the most recent analysis of Vie l et al. (2007), if 
o ne uses only the hig h-resolution high-redshift Lyman-a spectra 
of iBecker et alj ( |2007|) then one finds the lower bound on the ster- 
ile neutrino DM mass in the DW scenario to be ra, > 5.6 keV, 
which is in contradiction with our current upper bound m, < 4 keV 
(but would have left a narrow allowed window for in^ if one had 
used the previous bound m, < 6 keV o f I Asaka et al.ll2007l) . If one 
takes into account the low-resolution SDSS Lvman-O' dataset (Mc- 
Donald et al. 2006), used in Seljak et al. (2006); Viel et al. (200^), 
this contradiction becomes much stronger. Although the Lyman-a 
method relies on a very complicated analysis with (???) some un- 
known systematic uncertainties, it seems that the model in which 
all of the DM is produces through the DW scenario is ruled out. 

However, there is another way to produce the sterile neutrino 
through oscillations with active ne utrinos (resonant p roduction in 
the presence of lepton asvmmetries. lShi & FuUeil 19991 (SF)). In this 
case, one qualitatively expects that the results of the Lyman-o- anal- 
ysis can be lowered by a significant amount, as for the same mass, 
the mean velocity (free-streaming length) in the SF model can be 
much lower than in the DW model. However, as sterile neutrinos 
are produced in the non-equilibrium way and their spectrum dif- 
fers significantly from the thermal one, the actual Lyman-o- bounds 
may depend not only on the free-streaming but also on the detailed 
shape of the spectrum. The detailed analysis of the SF production 
and corresponding re-analysis of the Lyman-ff data is needed. Cur- 
rently, the SF mechanism is not ruled out. 

Finally, there is also the possibility of production of the sterile 
neutrino DM through the d ecav of the light inflaton (Shaposhnikov 
& Tkachev .2006) , which cannot be ruled out by X-ray observa- 
tions. 

Therefore, the sterile neutrino remains a viable and interesting 
DM candidate, which can be either warm or cold. One of the most 
interesting ranges of parameters is that of low masses, which is also 
in the pot ential reach o f laborato rv experiments ( Bezrukov & Sha- 
pos hnikovl2007h and w ill be pr obed with future X -ray spectrome- 
ters ("Bova rskv et al.l2006ai ; lden Herder et alj2007l) .^' However, the 
search for the sterile neutrino DM signal in all energy ranges above 
Tremaine-Gunn limit should also be conducted. 
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